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Recently, Yu, Brown, and Chuang [Phys. Rev. A 71, 032341 (2005)] investigated the entangle- 
ment attainable from unitary transformed thermal states in liquid-state nuclear magnetic resonance 
(NMR). Their research gave an insight into the role of the entanglement in a liquid-state NMR quan- 
tum computer. Moreover, they attempted to reveal the role of mixed-state entanglement in quantum 
computing. However, they assumed that the Zeeman energy of each nuclear spin which corresponds 
to a qubit takes a common value for all; there is no chemical shift. In this paper, we research a 
model with the chemical shifts and analytically derive the physical parameter region where unitary 
transformed thermal states are entangled, by the positive partial transposition (PPT) criterion with 
respect to any bipartition. We examine the effect of the chemical shifts on the boundary between 
the separability and the nonseparability, and find it is negligible. 

PACS numbers: 03.67.Mn, 03.67.Lx 



I. INTRODUCTION 

Quantum mechanics has very different conceptual and 
mathematical characters from classical mechanics. The 
superposition and the entanglement (nonseparability) of 
quantum states make the difference clear. Recently, these 
characterizations have been utilized for a newly develop- 
ing quantum technology. Actually, quantum entangle- 
ment is deeply related to quantum information process- 
ing (QIP) Q. 

The role of entanglement in quantum computing has 
been researched from various viewpoints. In particular, 
the generation of entangled states by quantum dynam- 
ics, or quantum algorithms, will give us useful infor- 
mation on the above task d 9. Recently, Yu, Br own, 
and Chuang [3| investigated the entanglement of uni- 
tary transformed thermal states in a liquid-state nuclear 
magnetic resonance (NMR) quantum computer. Such 
states are defined as the density matrices transformed 
from the thermal state in it by a specific class of uni- 
tary operators. We call the unitary transformations 
Bell-transformations. Their definition will be explained 
in Sec. Ill Bl Furthermore, we call such states Bell- 
transformed thermal states. The thermal state in Ref. Q 
is separable (i.e., there is no quantum correlation) Q and 
characterized by two physical parameters: The one is the 
number of qubits, and the other is a measure of the state's 
polarization. The authors in Ref. @ studied two kinds 
of Bell-transformations. One of their central interests 
is the difference between the Bell-transformed thermal 
states and effective pure states. The effective pure state 
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is the convex sum of the identity operator and a pure 
state, and a typical one is used in the current liquid 
state NMR quantum computer 0, 0| • They concluded 
that the former should be more easily entangled than the 
latter; the Bell-transformed thermal states can be entan- 
gled even in the physical parameter region where effective 
pure states are separable. 

Their research is very important for the following three 
reasons. First, it gives an insight into the role of the en- 
tanglement in a liquid-state NMR quantum computer. 
Braunstein et al. [8] pointed out that the effective pure 
states should be almost separable in the current liquid 
state NMR experiments. After their study, various stud- 
ies on the role of entanglement in liquid-state NMR 
quantum computing were done (see Ref. [3[ for additional 
references) . One should note that the most natural quan- 
tum state is the thermal one in the liquid-state NMR 
and an effective pure state is constructed from it by a 
sequential operation of several quantum gates. There- 
fore, the authors in Ref. 0] investigated the entanglement 
of a more elementary state than effective pure states in 
liquid-state NMR quantum computers. Secondly, they 
attempted to reveal the role of mixed-state entanglement 
in quantum computing. Its evaluation and the mean- 
ing of mixed states will be more subtle than the case 
of pure states. Nevertheless, they can be characterized 
from the viewpoint of quantum communication. In par- 
ticular, distillability is important [t| [l(| EH- If a mixed 
state is distillable, it is useful for QIP (e.g., quantum 
teleportation); we can distill a maximal entangled state 
by using a number of the copies, and local operations 
and classical communication (LOCC). Accordingly, their 
research could lead to the alternative understanding of 
mixed-state entanglement from the viewpoint of quan- 
tum computing. Finally, as mentioned in the beginning 
of this section, their work is related to the generation 
of entanglement by quantum dynamics. Actually, they 
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found the difference of the entanglement generation be- 
tween the two Bell-transformations. 

In this paper, we analytically derive the physical 
parameter region where the Bell-transformed thermal 
states are entangled in a more general case; specifically, 
the effect of the chemical shift [3, [13 is included in our 
Hamiltonian. The authors in Ref. [3] assumed that the 
Zeeman energy of each nuclear spin which corresponds 
to a qubit is a common value for all. However, it implies 
that one can't access the individual qubit; it is not a re- 
alistic model. Including such an effect will be necessary 
to evaluate the attainability of entanglement from the 
thermal state in the experiments. Moreover, the analysis 
taking account of the chemical shift will reveal how the 
difference between the Bell-transformations reflect on the 
physical parameter region where the Bell-transformed 
thermal states are entangled. Our method of evaluat- 
ing the entanglement is simple and straightforward: Pos- 
itive partial transposition (PPT) criterion [13J, [lj] . We 
analytically calculate the eigenvalues of the partial trans- 
posed Bell-transformed thermal states and find the least 
ones. Here, we would like to emphasize that, in Ref. 
the Diir-Cirac classification [15| has been applied to the 
evaluation of the entanglement about two types of Bell- 
transformed thermal states, but it doesn't necessarily 
work in all cases. 

This paper is organized as follows. In SecUH we in- 
troduce the Hamiltonian and the Bell-transformations. 
As in Ref. [3|], we concentrate on two unitary operators 
of Bell-transformations. In Scc. lIIIl in the first, we ex- 
plain how to specify an individual bipartition and briefly 
review the PPT criterion. Next, we explain the Bell- 
diagonal state [H, H3] in iV-qubit system, which plays a 
central role in this paper. Then, we show the main results 
in Sec. lIII Dl The sufficient conditions for the nonsepara- 
bility of the Bell-transformed thermal states with respect 
to any bipartition. In Sec. lIVl we show the necessary con- 
ditions for the full separability and the full distillability. 
In particular, we show the necessary and sufficient condi- 
tions for their full distillability when there is no chemical 
shift. In Sec. El we examine the effect of the chemical 
shift on the boundary between the separability and the 
nonseparability determined by the PPT criterion. We 
summarize our results in Sec. lVIl Furthermore, we briefly 
review the Diir-Cirac classification and show some exam- 
ples which this method doesn't work in Appendix. [Al 



II. MODEL 

A. Hamiltonian 

In a liquid-state NMR quantum computer, the qubit 
is the nuclear spin in the molecule. We assume the num- 
ber of the qubits is N in one molecule. The dipolc- 
dipole interaction between the molecules in solutions 
are negligible because they randomly collide with each 
other [7[. Therefore, we concentrate on the internal de- 



grees of freedom (i.e., the nuclear spin) of one molecule. 
Let us define the computational bases as |0)i and 
(i{0\0)i = 1, = I, and ,(011), = 0). The 

subscription i(= 1,2...N) is the label of the qubit. 
Furthermore, we introduce the following standard op- 
erators: I t = |Q),<0| + Z t = |0)i(0| - 

X t = \0)i(l\ + \l)i(0\, and Hi = (Z t + X t )/y/2. We 
analyze the Hamiltonian 

N h 

i=l 

where hvi is the Zeeman energy of the zth qubit. Note 
that we neglect the J-coupling between the neighbor- 
ing qubits because it is much smaller than the Zeeman 
energy in a liquid-state NMR 0, [Hj]. In this paper, we 
don't mention the relationship between the above compu- 
tational basis and the physical spin state, and the Hamil- 
tonian (fTJ) is a mathematical model. We assume hvi is 
positive when we derive our main results. However, the 
model includes the important physical effect; every value 
of hvi is different from each other, due to the chemical 
shift. The difference allows us to access each qubit indi- 
vidually. In Ref. [3J, every hvi is a common value. 

B. Bell— transformed thermal states 

Let us consider the separable state which is character- 
ized by a set of physical parameters. Our interest is the 
parameter region in which the state transformed from a 
separable one by quantum gates is a nonseparable one. 
Therefore, we have to specify a suitable initial separable 
state and quantum gates which generate the entangle- 
ment in a liquid-state NMR quantum computer. 

First, we refer to the initial separable state. The sys- 
tem in a liquid-state NMR experiment is a thermal equi- 
librium state with high temperature. Accordingly, the 
most natural choice for the initial state is the thermal 
state 



In our model, the Hamiltonian TL is given by Eq. |T]). 
Here, (3 and Z — tre - ^ are the inverse temperature and 
the partition function, respectively. Notice that state $2$ 
is a separable state with respect to any bipartition of 
the system. The physical parameters in Eq. © are the 
number of the qubit N, the Zeeman energy of the ith 
qubit hvi, and the inverse temperature /3, which are re- 
garded as the free parameters. Actually, the values of the 
parameters are restricted by the experiments for a liquid 
state NMR quantum computing. For the comparison 
with Ref. [H, we introduce the parameter a; = (3hvij2, 
which is called a measure of the state's polarization in 
the reference. Hereafter, we simply call it polarization. 
The physical meaning of a, is quite clear. Let us con- 
sider the expectation value of the z component of the 
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total spin operator J z (= X^=i Zi/2) with respect to p t h- 
m z = tr(J z pti,) = — ^^ 1 (tanhQj)/2. We assume the 
value of N is fixed. The value of \m z \ becomes larger as 
cti increases; the system is strongly polarized in the direc- 
tion of the z axis. Summarizing the above argument, we 
can say that the initial thermal state (J2j) is characterized 
by the number of qubits N and the polarization ai. 

Next, we explain the unitary operators for generating 
entanglement. We will have to investigate all types of 
the quantum gates which are considered as the generator 
of entanglement and essential parts in a quantum algo- 
rithm. However, this task will be very difficult. In this 
paper, as in Ref. we concentrate on the following two 
unitary operators: the controlled-NOT-Hadamard (CH) 
transformation, Uqh = Uf an (Hi®I), and the CH-fanout 
transformation, Ucf = UcnU{ an . The fanout gate J7f an 
is defined as follows: U lan = |0)i(0| <8> I + ® X, 

where I = ® i=2 hi X = ®^1 2 X- The above uni- 
tary operators are examples in quantum gates, but they 
include the essential quantum gates for the generation 
of entanglement: controlled-NOT gates and Hadamard 
gates. In this case, they generate the entanglement be- 
tween the first qubit and the remaining qubits. When 
we consider, for example, the case of N = 3 and that 
the initial state is |000) = |0)i <8> |0) 2 ® |0) 3 , we obtain 
the following results: C/ C h|000) = (|000) + |lll))/v / 2 and 
C/ C f|000) = (|000) + |111))/V2. 

In summary, we examine the entanglement of the states 

PCH = f/cHjOth^CH' ( 3 ) 

Pcf = U C FPthUl F . (4) 

We call the above states the Bell-transformed thermal 
states. The above states are called the Bell-transformed 
thermal states. The initial state pth is separable with 
any bipartition of the system for the arbitrary values of 
N and {ai}g v 



III. ENTANGLEMENT OF 
BELL-TRANSFORMED THERMAL STATES 
WITH CHEMICAL SHIFT 

A. Specification of a bipartition 

In order to study the entanglement of a system, it 
is necessary to specify the way to divide it into two 
parts. We divide the A-qubit system into two subsys- 
tems, party A and party B, in the following 3]. First, let 
us consider a set of binary numbers, {ki}f =1 (hi = 0, 1). 
When fcj = 0, let the zth qubit be in party A. On the other 
hand, when fcj = 1, it is in party B. We always set k\ — 0; 
the first qubit is always in party A. For convenience, we 
introduce an integer k = ^2^ =2 ^i^ ~ 2 ■ Therefore, a par- 
tition is specified if an integer k(G [1, 2 N ~~ 1 — 1]) is chosen; 
we call such a partition the bipartition k. Let us choose, 
for instance, k = 4 in the case of N = 4 (i.e., fc 2 = 0, 



ks = 0, and k^ = 1). The elements of party A are the 
1st, 2nd, and 3rd qubit, and party B contains only the 
4th qubit. 

B. PPT criterion 

The PPT criterion is the simple and computable way 
to investigate entanglement. We briefly recapitulate it. 
Let us consider a density matrix p in a quantum system 
with finite dimension d. The total system is divided into 
two subsystems, system A and system B. Introducing an 
orthonormal basis of the system A, {|i/j)A}f=i and the 
system B, {Iv/Ob}^! j we can expand the density matrix 
p as follows: 

P= ^2 C(ik\jl) \ui) A (uj\ ® K)b(«*|> (5) 

i,j=l k,l=l 

where C(ik\jl) is a complex number and d — dj^dB ■ Next, 
using Eq. ([5|) , we define the partial transposition of the 
density matrix with respect to the system B as 

p Tb =J2 E c(iZ|jAOk>A<^|®k>BM (6) 

i,j=l k, 1=1 

Then, we calculate the eigenvalues of p Tu and investigate 
their positivity. If all eigenvalues of Eq. ([6]) are positive 
(i.e., p TB > 0), it is called a density matrix with PPT. 
On the other hand, if at least one eigenvalue of it is neg- 
ative, it is called a density matrix with negative partial 
transposition (NPT). The most important thing is that 
there is the following criterion (i.e., PPT criterion): 

p: separable => p: PPT, (7) 

or, equivalently 

p: NPT -> p: entangled (nonseparable). (8) 
Moreover, the following statement is also known [ToL [TlT|: 
p: distillable => p: NPT. (9) 

C. Bell— diagonal states 

Before showing the our results, we explain the special 
class of a density matrix, the Bell-diagonal state [l6|, [ijj • 
It plays a central role in later discussion. 

To begin, let us explain thegeneralized Greenberger- 
Horne-Zeilinger (GHZ) state [J OS m tne AT-qubit sys- 
tem 

|^) = ^(|0i)±|lj)) (0<j<2 w - 1 -l), (10) 

where j = YliLz ii2 l ~ 2 f° r the binary number (= 0, 1), 
|0j) - |0)i ® ®^ 2 \j t )i and = |l>i O ®t 2 |1 - 
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The symbol j means a bit-flip of j: j = 2 JV_1 — 1 — j. 
We can easily find the generalized GHZ states are the 
elements of an orthonormal basis of the Hilbert space 
corresponding to the A-qubit system. 

We introduce the following density matrix: 

2 «-l_ 1 

pbd= J2 Kl*i>W+ w i1*7><*7l)' C 11 ) 

where uf = (^flpBB^f) and E^" 1 H + U i) = 
1. Equation pip is a Bell-diagonal state in an A-qubit 
system. We will show that pcH and pcf take the form of 
Eq. pip in the next subsection. 

The reason why we introduce the Bell-diagonal state 
is that we can easily obtain its partial transposition [13] ■ 
We confirm this in the following procedure. First, 
let us consider a bipartition k. Secondly, we repre- 
sent \^f)(^f\ in the computational basis: |^)(^| = 
(|0i)(0j| ± |0j)<T| ± |H)(0j| + |lj)(lj|)/2. The diag- 
onal parts of |*J)(*J| are |0j)(0j| and and 
are invariant under the partial transposition with re- 
spect to party B. The off-diagonal ones are |0j)(lj| and 
|lj)(0j| because j ^ j. If the ith qubit is in party B 
(i.e., ki = 1), the binary number j of the off-diagonal 
parts is transformed into ji + 1(= ji + ki) modulo 2 
(e.g., 0+1 = 1 and 1 + 1 = 0) by partial transpo- 
sition with respect to party B. On the other hand, if 
the ith qubit is in the party A (i.e., ki = 0), the cor- 
responding ji is unchanged; ji = ji + fcj. As a result, 
we obtain the following expression: (|\l/^ = )(^'j b |) TB = 

><*^*M*7 efc K*7 efc l)/2, 
f _ 2 «»* \-i=ji + h mod 2). Accord- 
ingly, the Bell-diagonal state partially transposed with 
respect to party B is given by 



wim\+\*jm\±\* i&t 

where j®* = E^a z «2 i_2 (U 



2 n-i_ 1 

pbd = E (^ + l*/><*;i+M7l*7X*7l)'( 12 ) 

3=0 



where 



± 



(13) 



where < j, j' < 2 Ar ~ 1 — 1. Consequently, both /9ch and 
Pcf are Bell-diagonal states. 

Let us define the mean values of the polarization of 
party A, £ and party B, rj for given j and k as follows: 



* = AT" E ^ 
N — w *—f 



oti, r\ 



= I^(-l)iia i , (17) 



ieB k 



where A k = {i G Z; ki = 0, 1 < i < N}, B k = {i e 
= 1, 1 < i < A}. In Ea. ([T7]). we conventionally 
assign zero to j\. The number w is the total number of 
the elements of party B (1 < w < N— 1). In other words, 
it is the hamming weight of k (i.e., the number of one in 
{h}iLi)- Each of £ and r] is a function of j if A, {ai}^L 1: 
and k are fixed. 

Then, we calculate the eigenvalues of and p^p 
with respect to the bipartition k. Now that we know 
the general expression (flU)) for the partially transposed 
Bell-diagonal state, we easily obtain the desired results. 
For a given bipartition k, the eigenvalue of is given 
by 



(18) 



where 

n% H (r)) = e- wr > + tanhai e wr >. (19) 
Similary, the eigenvalue of p^,p is given by 

1 



where 



MCF, j = ~f n C F (£> 



n± F (£, r?) = e^-^cosh^) 



(20) 



±e ±(Jv-«)« sinh ( W7? ). (2i) 

In general, the relative signs between c^s can be differ- 
ent. For the latter discussion, we evaluate the range of £ 
and rj for the given {cti}f =l and k. We rewrite Eq. (JTTJ) 
as 

a = E (-^"M* = - E (- 1 ) i<+ " 4 i«*i. 

/V — in — ' in — ' 



ieA k 



iEB k 



D. Sufficient conditions for nonseparability 

First of all, let us show Eqs. ([3|) and |4]) are just the 
Bell-diagonal state. Using the standard relations (Hi ® 
Wi) = (\0j) + M)/V2, (H^I)\lj) = (|0j)-|lj))/V2, 
^fan|0j) = \0j), and C/f an |lj) = we obtain the fol- 
lowing results: 



(*f|p CH |*£) 



^ZZl e T"i e -E™ 2 (-i) 31 ^ 



j irvr , r , z 

± I IVT>T\ /iTr± 



^ZZl e T"i e =FE™ 2 (-i) 31 ^ 



(*>ch|*J) = (*f |pcf|^) = 0, 



(14) 

(15) 
(16) 



where Sj = for positive eti and s» = 1 for negative a» 
We readily obtain the inequalities 

/V — 9/ ^ ^ 7/T ' 



The condition for £ = £*(= £ m aa;) is easily found: ji+Si = 
mod 2 for any i € A^. Similarly, the conditions for 77 = 
»?*(= W) and r] = -r] 1f (= Tj mm ) are given by jj + Sj = 
mod 2 and j'j+Sj = 1 mod 2 for any i £ B k , respectively. 
Notice that the minimum value of £ is not always — 
because ji = 0; it is -£* + [1 + (-l) Sl ]|ai|/(iV - io) = 

£>rnin ■ 
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Now, we analytically derive the sufficient conditions for 
the nonseparability of pen and pep- What is needed is 
that we search for the minimum values of the eigenvalues 
which can be negative. 

First, we examine Eq. (|18p . In order to discuss defi- 
nitely, we assume u\ is positive for a while. We find that 
the eigenvalue p^, H . is always positive because all fac- 
tors of the right hand side are positive. Then we focus 
on /zJ H .. The positivity is determined by the value of 
n CH (v) ■ Notice that it is a monotonic decreasing function 
of 77, because dn^^rj) / drj < 0. Therefore, the minimum 
value of n£ H (7?) is given by nJ H (rj max ) = 71^(77*) = 
e -«"7* — tanhai e wr> * . We examine the case of negative 
oti in turn. In the case, the value of 71^(77) is impor- 
tant for the examination of the positivity of pen- We 
can readily check that the minimum value of 77^(77) is 
given by n^ H (r] max ) = n CH (r;») = e~ W11 ' -tanh |ai| e wr i*. 
Summarizing the above argument, we can say that pch 
is NPT with respect to the bipartition k if and only if 

n£ H (>7*)<0 e- 2w, i* < tanh|ai|. (23) 

Secondly, let us consider Eq. f20]). We concentrate 
on the behavior of n^ F (^, 77) because Z > 0. First, 
we investigate nJ F (£, 77). Notice that the value is 
always positive for 7/ > or £, 77 < 0. Here- 
after, we consider the case of £ > and 77 < 
0. In this case, we find <9nJ F (£, 77) /<9£ < and 
dn^, F (^, rf)/dr] > 0. Therefore, the minimum value of 
n CF(£> V) is given by n^, F (£ max , r\ min ) = n£ F (£*, -7?*) = 
e -(iV-iD)£« cosh(7i;?7») — e( N ~ w ^* sinh(u>?7»). Next, we 
consider Uq f (£, 77). It should be noted that the fol- 
lowing relation is fulfilled: 71q F (£, 77) = n^ F (— £, —77). 
Therefore, we readily obtain the information on 
tt-C! F (£, v)- Through the above arguments, the mini- 
mum value of 7ip F (£, 77) is given by n^ F (^ min , r\ max ) = 

ncAZmin, V*) = e-W-vX'+W-WM coah(wri.) - 
e (iV— u,)£»-[i+(-i)-i]|« a | sin h( w ^). On the other hand, 

n CF r 7*) i s clearly greater than or equal to 

7t.q F (£*, —?7*)- Consequently, pcf is NPT with respect 
to the bipartition k if and only if 

77-(t F (£*, -?7*) < 

cosh [(N - w)£* - wrj*] < sinh (Not) , (24) 

where iVa = Eili l Q d = (N — + w^*- 

We summarize the sufficient conditions for the nonsep- 
arability of pcH and pcF with respect to the bipartition 
k in TableU Note that the information of the bipartition 
k is included in and 77, through Ak and Bk- The defi- 
nitions of Ak and Bk are explained below Eq. (fTT|) . With 
respect to the bipartition fc, pch is an entangled state 
if the inequality (f2"3"|) is fulfilled, and pcF is an entangled 
one if the inequality ([24]) is fulfilled. The sufficient condi- 
tion for the nonseparability of pch is given by ol\ and the 
mean value of the magnitude of polarization for party B, 
77*; it doesn't depend on on (i £ Ak) in party A, except 



TABLE I: Necessary and sufficient conditions for NPT with 
respect to the bipartition k for each Bell-transformation, 
where = J2 isAk \<*i\/( N ~ w )> V* = E ieSfc \<Xi\/w, and 



Bell-transformation 


condition for NPT 


CH 


e -2w V , < tanh | ai | 


CH-fanout 


cosh [(N — to)£* — wq*\ < sinh (Na) 



for ol\. On the other hand, the one of pcF is determined 
by £* and 77* ; it depends on the mean values of the mag- 
nitude of polarization for party A and party B. 

It should be noticed that if every on is a common value, 
we can readily check that both sufficient conditions (|2"3"|) 
and (124p are equivalent to the corresponding results in 
Ref. % 

E. Characterization of the Bell— transformations 

It is important for the deep understanding of quan- 
tum algorithm to characterize the property of quantum 
dynamics in terms of entanglement. Let us consider 
the ability and mechanism of the Bell-transformations 
C/ch and Uqf to generate entanglement from the ther- 
mal state. 

The authors in Ref. [|| discussed this ability without 
taking account of the chemical shift. They concluded 
that Ucf is a more effective Bell-transformation than 
C/cHj because the parameter region, in which pqf is en- 
tangled with respect to a given bipartition, is wider than 
the corresponding one of pen- Such a result is quite nat- 
ural because the number of the controlled-NOT gates in 
Ucf is twice as many as in E/ch- 

In Sec. lIIIDl the analysis of the sufficient conditions 
for the nonseparability of the Bell-transformed thermal 
states with the chemical shift reveal the difference be- 
tween Uch and Ucf- We assume that the maximum 
value of 77, 77* is given with respect to a bipartition k. 
In the case, let us consider the necessary information 
of party A to examine the entanglement of the Bell- 
transformed thermal states by the PPT criterion. As 
for pch, from Eq. (|23[) . we find that the value of a% is 
only needed; the local information of party A is required. 
On the other hand, we have to know the maximum value 
of £, to investigate whether pcF is entangled or not, 
according to Eq. (|2"4")k the global information of party A 
is required. These observations can lead to the under- 
standing how the Bell-transformations make party A and 
party B entangled. 

IV. FULL SEPARABILITY AND FULL 
DISTILLABILITY 

Let us consider a ./V-particle system. A state of this 
system is called fully separable (or, TV-separable) if the 
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corresponding density matrix p can be written as a con- 
vex combination of direct product states: 

i 3=1 V i ) 

(i) 

where p~ is the density matrix on t he p artial Hilbert 
space corresponding to the jth partile [Hf, [l8| . One can 
easily check that a density matrix has PPT with respect 
to any bipartition of the system if it is fully separable. 
On the other hand, we call p fully distillable if it is distil- 
lable with respect to any bipartition (i.e., we can create 
a maximal entangled pair between qubits in party A and 
party B by a number of the copies and LOCC), accord- 
ing to Ref. [3j . Through statement ^ , we can readily 
show that a density matrix has NPT with respect to any 
bipartition, if it is fully distillable. 

As in Ref. [|, we examine the full separability and the 
full distillability of pcH and Pcf- We summarize the 
necessary condition for the separability of pcH and pcf- 
According to Eqs. (j2"3")l and (|2"4")l . the necessary conditions 
for the separability of pcH and pcf, respectively, with 
respect to the bipartition k are given by 

e - 2wT >> > tanh|ai|, (26) 
cosh [(N - w% - > sinh (No) . (27) 

Moreover, from Eq. ([9]), we find that the necessary con- 
ditions for the distillability of pen and pcf with respect 
to the bipartition k are given by Eqs. (j2"3"|) and (f2"4"|) , re- 
spectively. Note that the equalities of Eqs. (f2"6"| and (j2"T|) 
give the boundaries between the separability and the dis- 
tillability. 

Let us consider the full separability and the full 
distillability of pen- The necessary condition for 
the full separability and the full distillability is given 
by mm k (e- 2w ^) > tanh | ail and max fc (e~ 2wr >*) < 
tanh |ai|, respectively. In conclusion, we obtain the nec- 
essary conditions for the full separability and the full dis- 
tillability of pea as follows: 

pen- fully separable e ~ 2bmax > tanh |ai|, (28) 
pen- fully distillable => e ~ 2bmin < tanh|ai|, (29) 

where b max = max. k (wn*) = Na - \ai\ and b mm = 
minfe(w77*) = minj^i(|oj |). The value of 1/777* takes 
the maximum value b max when k = 2 N ~ 1 — 1 (i.e., 
W = N — 1). On the other hand, it takes the minimum 
value brain when ki corresponding to the minimum value 
of ajS is 1 and the remainders are (i.e., w = 1). 
Similarly, we obtain the following results for Pcf'- 

Pcf ■ fully separable cosh d m in > sinh (Na) , (30) 
Pcf- fully distillable coshd moa! < sinh(Aa) , (31) 

where d max = max fc |(A - w)^ - wr)*\ and d min = 
minfc I (N — w)£* — wr]* | . The value of | ( N — w)^ — wq* | 
takes the maximum value, d ma x, when the difference be- 
tween (N — w)t;* and wrj* is the largest. When a, = 



a(> 0) for any i, the condition for d max is quite simple. 
The value of d max is (N — 2)a, where k = 2 7V_1 — 1 (i.e., 
w = N— 1) or k = 1 (i.e., w = N—l). On the other hand, 
it takes the minimum value, d m i n when (N — w)^ is the 
closest value to wry*. If every ai is a common value, the 
value of dmin is 0, where w = [N/2\ . Here, the symbol 
[a;] means the greatest integer that is less than or equal 
to x G E. 

Our results (|28|) - ((3"Tj) are the necessary conditions 
for the full separability or the full distillability. When 
a t = a(> 0) for any i, the authors in Ref. [|[ showed the 
sufficient condition for the full distillability, through the 
statements proved in Ref. [15j . In this case, combinating 
our results (f2"9"|) and (|3T|) with theirs, we can obtain the 
following important results: 

Pch: fully distillable <=4> e~~ 2a < tanha, (32) 
Pcf: fully distillable 4=^ tanha > e" 2(Ar " 1)a , (33) 

(We can find that the necessary condition for the full 
distillability of pcf is cosh[(iV — 2)a] < sinh(TVa) from 
Eq. (|31[) and this inequality is equal to the corresponding 
expression in Eq. (j33|) after a short calculation). There- 
fore, we obtain the complete physical parameter regions 
in which pch and pcf can be useful for QIP. 

V. BOUNDARY BETWEEN SEPARABILITY 
AND NONSEPARABILITY 

The equalities in the inequalities (f2"6")l and ([2T| imply 
the boundary between the separability and the nonsepa- 
rability in term of the PPT criterion. We investigate the 
effect of the chemical shift on such boundaries. Hereafter, 
we assume all aiS are positive, for the sake of simplicity. 

Through the above argument, we easily obtain the fol- 
lowing expression of the boundary for pen ■ 

e - 2wr >* =tanho!i. (34) 

Similarly, the boundary for pcf is given by 

cosh [(N — w)£* — wr]*} = sinh (Na) . (35) 

To compare the boundaries with the chemical shift to 
those without it, we try a toy model for a%. The po- 
larization ai divides into a reference value a(> 0) and 
the deviation 5ai from it: at — a + 8ai. Here, we re- 
gard Xi = Sail a as a uniform random variable in [—5, 5] 
(0 < S < 1). Therefore, in our toy model, the polariza- 
tion of the ith qubit is given by «i = a(l + s t ). Fur- 
thermore, the value of a is regarded as the mean value 
of polarization. 

Let us explain how to calculate the boundaries. In the 
first, we choose the value of S; actually, S = 0.1, 0.01, 
and 0. Then, a sequence of random numbers in [—6, S] is 
generated by the Mersenne Twister [l|. Next, we specify 
a bipartition. For the sake of simplicity, we assume the 
first N — w qubits are in the party A and the latter w 
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qubits are in party B. Finally, we numerically calculate 
the value of a satisfying with Eqs. (j3"4"|) or (|35[) for a given 
N by the bisection method. 
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FIG. 1: Boundary between the separability and the nonsep- 
arability in term of the PPT criterion. The solid line is for 
8 — 0.1, the broken line is for 8 = 0.01, and the broken 
dotted line is for 8 = 0; actually, the difference between 
these is very small (see the insets), (a) The CH transfor- 
mation (k = 2 iv ~ 1 — 1). (b) The CH-fanout transformation 
(k = 2 N - 1 -2^- 1 ). 

In Figs.QJa) and (b), we show the boundaries on 
(a, N) plan for a specific sequence of X{. The horizontal 
axis is the common logarithm of the inverse of a. The 
larger value of log 10 or 1 corresponds to the case of the 
higher temperature. The longitudinal axis is the number 
of qubits. The solid line, the broken line, and the bro- 
ken dotted line correspond to the boundaries for 5 = 0.1, 
8 = 0.01, and 8 — 0, respectively. We show the results 
for pen with respect to a bipartition k = 2 N ~ 1 — 1 (i.e., 
w = N — 1) in Fig.QJa), and the case of pcf with respect 
to a bipartition k = 2 N ~ 1 -2^ N / 2 i~ 1 (i.e., w = |_iV/2j) in 
Fig.QJb). The left sides of those lines are the parameter 
regions where pch and pcf are nonseparable with respect 
to the corresponding bipartition. We find that the clear 
distinction among the different values of 8 is invisible in 
both Figs. [I] (a) and (b); for example, denoting the value 
of a on the boundary as ab(5) for a given 8, we can find 
that |log 10 [a b (0)/a b (S)} \ < 10~ 2 for 8 = 0.1 and 0.01, 
as N is large. In addition, we can find a similar behav- 
ior even if we change the sequence of random variables 
and the kind of bipartitions. Consequently, the effect of 
the chemical shift on the boundary between the separa- 
bility and the nonseparability is negligible in our model 



for {cti}f =1 . This result implies that one have only to 
examine the number of qubits N and the mean value of 
polarization a for the determination of the entanglement 
of the Bell-transformed thermal states (J3]) and (H|) . 

Let us consider the reason why the effect of the chem- 
ical shift is negligible. The boundaries for both pch and 
Pcf between the separability and the nonseparability are 
given by Eqs. (|3"4"|) and respectively. Those equa- 

tions are mainly determined by the mean values of the 
polarization of party A, and party B, 77*. Therefore, 
the deviation from the mean value of o^s is not impor- 
tant for the determination of the boundaries. However, 
the model for a; is quite simple; the distribution of a; 
is uniform and random. We have also assumed that ajS 
are positive in this section. Therefore, it is necessary to 
investigate a more general and realistic model for a,. 



VI. SUMMARY 

We have analytically derived the physical parameter 
region where the Bell-transformed thermal states are 
entangled in the presence of the chemical shift, by the 
use of the PPT criterion with respect to any bipartition. 
Two kinds of Bell-transformations, the CH transforma- 
tion Uca and the CH-fanout transformation £/cf, have 
been examined, as in Ref. With respect to the bipar- 
tition k, pen is an entangled state if the inequality (|23p 
is fulfilled, and pcf is an entangled one if the inequal- 
ity ([24]) is fulfilled. We summarize our results in Table 
HJ If the every a; is a common value, our results are 
equal to the corresponding ones in Ref. There exists 
an obvious differences between Eqs. (p2"3")l and with 
respect to their dependence on on. The sufficient con- 
dition for the nonseparability of pch is given by ot\ and 
the mean value of \on\ for party B, 77*; it doesn't depend 
on cti (i € Ak) in party A, except for a\. On the other 
hand, the sufficient condition of pcf is determined by 
and 77*, the mean values of \ai\ for party A and party B, 
respectively. 

Next, we have derived the necessary conditions for the 
full separability and the full distillability through the 
above analytical results. On the other hand, the authors 
in Ref. [|[ obtained the sufficient conditions when all of 
aiS are equal. We have obtained the necessary and suffi- 
cient conditions for the full distillability of pen and pcf 
in this case, combinating our results with theirs. Accord- 
ingly, we can completely characterize the physical param- 
eter regions in which those unitary transformed thermal 
states can be useful for QIP, if there is no chemical shift. 
When the chemical shift exists, we will have to examine 
the sufficiency of our results. 

Finally, we have investigated the effect of the chemical 
shift on the boundary between the separability and the 
nonseparability determined by the PPT criterion. We 
have shown such an effect on the boundaries should be 
negligible. Actually, one have only to know the number 
of qubits N and the mean value of the magnitude of po- 
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larization |aj| for the evaluation of the entanglement of 
the Bell-transformed thermal states ([3]) and (fj|. This 
result is quite natural, because the boundaries for both 
Pch and pcf between the separability and the nonsep- 
arability dominantly depend on the mean values of |oi| 
for party A, and party B, 77*. We will have to research 
a more general and realistic model for on for examining 
the effect of the chemical shifts on the boundary between 
the separability and the nonseparability. 

It is also important to examine the difference between 
Pch and pcf (ke., Uch and Ucf) in terms of the entan- 
glement generation. The authors in Ref. [I| concluded 
that Ucf is a more effective Bell-transformation than 
Uch because the parameter region in which pcf is fully 
distillable is wider than the corresponding one of pen- 
The analysis taking account of the chemical shift has re- 
vealed the difference between Uch and Ucf from another 
point of view. When the value of 77* is given with respect 
to a fixed bipartition k, the value of ot\ is necessary to 
examine the entanglement of pch ; the local information 
of party A is required. On the other hand, for Pcf, the 
mean value of \a.i\ for party A, is essential; the global 
information of party A is required. In this paper, we 
have obtained the analytical expressions of the eigenval- 
ues of Peg and pjp. Therefore, the results allow us to 
characterize Uch and Ucf in more detail; for instance, 
we can evaluate the negativity, which is an entanglement 
measure [13, [U| . We will show the results in the near 
future. 

Research on entanglement in liquid-state NMR in- 
volves various aspects of quantum information theory, for 
example, the role of mixed states in quantum comput- 
ing and the classification of entanglement. One should 
note that the achievable range of the physical parame- 
ters N and a^s is limited in a current liquid-state NMR 
experiment. Actually, it may be difficult to compare 
the theoretical results with the experiments. However, 
several experimental developments have been reported 
in liquid-state NMR, for example, the highly polarized 
initial states [22j and the number of qubits greater than 
ten [23j, l24[ . Furthermore, research on a solid-state NMR 
quantum computer [2^, [26j , which can relax the limita- 
tion of liquid-state NMR, has been developed steadily. 
Consequently, we expect that theoretical research on the 
entanglement in liquid-state NMR could be connected 
with these experiments in future. 
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APPENDIX A: NOTES ON DUR-CIRAC 
CLASSIFICATION 

The Dur-Cirac classification [l5| is a very effective way 
to evaluate the separability or the distillability of density 
matrices, of either pure or mixed states, in a multiqubit 
system, and it is widely used [1, [H, |27|. However, we 
have to take care of using it; we don't always obtain infor- 
mation on their entanglement from the method proposed 
in Ref. [15J- In the appendix, we show several examples in 
which the method doesn't work. Unfortunately, we don't 
know what kind of density matrices have the problem. 
Nevertheless, in the end of the appendix, we propose the 
prescription to solve this problem in specific examples. 

First of all, we summarize the strategy of the Dur- 
Cirac classification. In order to evaluate the entangle- 
ment of a quantum state by the PPT criterion, one must 
examine the positivity of the partial transposed density 
matrices. In general, such tasks will be difficult as the 
number of qubits becomes large and many choices of bi- 
partitions exist. On the other hand, if one uses the Dur- 
Cirac classification, it is only necessary to calculate some 
specific matrix elements of the state concerned. The main 
idea is that, using a sequence of local operations, one can 
transform an arbitrary density matrix p in a 7V-qubit 
system into the following state whose property of entan- 
glement is easily examined: 

PN = A+|*+)<#+| + Ao|*o)<* -| 

2 JV-i_! 

+ E x j(\*p(*i\+\*j)(*i\)- ( Ai ) 

The original density matrix p is related to pjv by the 
following expressions: A* = (^q |p|*o) anc ^ ^Aj = 
+ (Vj\p\yj). The property of p N is char- 
acterized as follows: 

p N : PPT with a bipartition k ^=> A < 2A fe , (A2) 

or 

p N : NPT with a bipartition k A > 2X k , (A3) 

where A = |Aq — A^~|. The most important key idea is 
that the entanglement does not increase through local 
operations. Accordingly, if p^ is a nonseparable state 
with a bipartition, p is also such a state. It should be 
noticed that we obtain information on the entanglement 
of p only if pj\r is a nonseparable state. 

Now, we show three examples. The first one is the case 
in which the Dur-Cirac classification works well. The 
remaining two are not such cases. Hereafter, we concen- 
trate on a two-qubit system. Therefore, the value of k 
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in Eqs. (|A2|) and (|A3[) is always 1. In the first, let us 
consider the following state: 



so = {l-f)\{h 



^)+/I^X*o + l, (A4) 



where —1/3 < / < 1. The above density matrix is called 
an isotropic state [ll[ . Directly using the PPT crite- 
rion [l3|, [3] > we readily find p iso is an entangled state 
if / > 1/3. On the other hand, we apply the Diir- 
Cirac method to Pi SD . We obtain the following results: 
(*+\p iso \V+) - (1 + 3/)/4, (¥„ \Pi.o\% ) = (1 - /)/4, 
and (*+|p iso |*+) + (^\p iso \^) = (1 - f)/2. Accord- 
ingly, using Eq. (jA3|) . we also find pi SO is an entangled 
one if / > 1/3. 

Next, we consider a slight different state from pi SO 



P,. 



(A5) 



where —1/3 < / < 1. Notice that the condition for 
the nonseparability of p' iso is the same one as pi SO ; p' iso 
is an entangled state if / > 1/3. This result is quite 
natural, because the state pi SO is transformed into p' iso 
by a local unitary operator; p' iso = (Zj ® X 2 )pi S o(Ii <8> 
X2)'. On the other hand, we obtain the eigenvalues 
in the form of Eq. (|A"l]) as follows: (*o IpLoI^o ) = 
(l-/)/4, (*o l/CI*o > = (l"/)/4, and (*+|^ |*+) + 
(*i l/4J*i) = (1 + /)/ 2 - Therefore, the value of A is 
less than 2Ai for the arbitrary value of / and the con- 
dition (|A2j) is satisfied. This might imply that p' iso is 
always separable. However, p' iso is local unitary equiva- 
lent to the entangled state p iso (f > 1/3). In conclusion, 
we can't obtain information on the entanglement of p' iso 
by the method in Ref. [HI], because the entangled state 
p' iso is transformed into a separable state by LOCO 

The third example is related to the task in this paper. 
We consider the case of N = 2 and otx = a? = ot in 
Eq. (J3]). We have known the condition for the nonsepara- 
bility of pcn'7 Pen is an entangled state if e~ 2a > tanha. 



However, the value of A is less than 2Ai for any a be- 
cause (*J|pch|*J) = e~^/Z, (* |p CH |*o) = V-Z 
and (*+|p CH |*^) + (*rlPc H |*r> = (1 + e 2a )/Z. Ac- 
cordingly, we can't also obtain information on the entan- 
glement of pen by the method in Ref. 15]. 



The authors in Ref. [|[ pointed out that the local op- 
erations can decrease the entanglement of the state con- 
cerned and an alternative method is necessary. In partic- 
ular, they discussed a more effective procedure for evalu- 
ating entanglement than the local operations in the Diir- 
Cirac classification in terms of the mathematical the- 
ory of majorization. Finally, we comment on a loop- 
hole in the Diir-Cirac classification from another point 
view: We attempt to construct a prescription to solve 
the problem. As has been mentioned, the isotropic state 
Pi SO is transformed into p' iso by the local unitary oper- 
ator Ii <g> X 2 . The method in Ref. works for the 
former but not for the latter. This suggests that, by 
suitable local unitary operators, the state for which the 
Diir-Cirac classification doesn't work, could be trans- 
formed into a proper one to which their method is appli- 
cable. Let us consider the above third example. When 
we use the local unitary operator I\ ® X 2 , the state con- 
cerned is transformed into p' CR — (Ii®X 2 )pcn(h®X 2 )^ , 
where (^\p'ck\*o) = 1/Z, CWchI*o > = ^l z ^ and 
(*t\p' CB \*t) + <*iVcHl*r> = (e- 2Q + l)/f Calcu- 
lating A and 2A 1; we can obtain the non-trivial expres- 
sion for p' CH ; for example, p' CH is NPT if and only if 
sinh(2a) > 1. If p' CH is entangled, the corresponding 
state pch is also considered to be entangled. Therefore, 
we obtain the sufficient condition of the nonseparability 
for pea, by the use of the Diir-Cirac classification; such 
a condition is just sinh(2a) > 1. This condition is equal 
to the one which is derived by the use of Eq. (|23|) . It 
is necessary to examine whether such a prescription is 
generalized or not. 
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